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ABSTRACT 
Some  sufficient  conditions  for  the  explicit  block  diagonal 
decomposition  of  block  matrices  are  developed.  These  applied  to 
matrices  corresponding  to  symmetric  and  regulcir  structures  of  finite 
size.  A  special  numerical  procedure  is  proposed  for  solving  linear 
systems  with  quasi  block  Toeplitz  matrix.  An  explicit  formula  for  the 
natural  frequencies  of  a  clamped  rectangular  plate  is  derived. 
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1 .   INTRODUCTION 

1.  Large  scale  matrices  (we  consider  them  as  block  matrices  A») 
usually  correspond  to  physical  structures  (models)  which  contain 
identical  or  periodically  repeated  substructures  (submodels),  for 
example,  matrices  corresponding  to  symmetric  or  regular  structures 
(models) .  In  some  cases  these  matrices  can  be  explicitly  block 
diagonalized  or  even  diagonalized.  They  are  the  subject  of  this  paper. 

2.  Introduce  some  definitions.  Let  A  be  nxn  non-defective  matrix 
(i.e.,  it  has  a  full  set  of  eigenvectors) 


A=  U,^]l^.,  (1.1) 


The  following  factorization  of  A  is  called  the  spectral  decomposition 


A  =  U  A  U"^  (1 .2) 


Here  U  is  the  modal  matrix  (its  columns  are  eigenvectors  of  A)  and  A  is 
the  spectral  matrix  (diagonal  matrix  with  eigenvalues  of  A), 

U=  [uij]?,j=i  .      A=  [Aj  JJ^T  (1.3) 

Special  brackets  [...J  and  one  subscript  are  used  for  diagonal  and  block 
diagonal  matrices:  f  X^  JV_i  -  [s.-^A^]^  ._.  ,  6^^  is  the  Kronecker 
delta.  We  shall  call  (1.2)  the  explicit  spectral  decomposition  of  A 
(1.1)  if  X^  and  u^^^  are  given  as  explicit  functions  in  a.,  and  n.  Let 
k^   be  a  block  matrix  of  order  mn  with  mxm  blocks 

A.  =  hij^.y-:  -   [[-ijaJa,.=ili,j=i  (1-^) 
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and  A»  be  another  block  matrix  with  the  same  elements  a^.^^^  grouped  in 
nxn  blocks 

~^*-lKX,.-^   =  [[^jox]?.j=i];..=i  (1.5) 

Both  matrices  are  similar 


A^^  =  pIa»P,  ,  (1.6) 


where  P^^  is  the  following  permutation  matrix  with  rectangular  mxn 
blocks 


^*  =  rf^ij^T];:J=i  ]?:?=!  (lt) 


Obviously 


p;^  =  pi  =  [[6ij^j^:M?:?=i  (1.8) 


The  following  block  diagonal  decomposition  of  A»  (1.^) 


A^  =  U»  A»  u;1  ,  (1.9) 


where  A,^  is  block  diagonal,   is  called  the  explicit  block  diagonal 
decomposition  of  A«  if  blocks  of  ^^   and  U,^  are  explicit  functions  in 
eigenvalues  and  eigenvectors  of  Aj;.  and  m,n. 
3.  Great  advantage  in  solving 

(A»-AI)x»  =  b»  (1.10) 


can  be  achieved  If  A«  possesses  the  explicit  block  diagonal 
decomposition.  The  necessary  condition  for  the  explicit  block  diagonal 
decomposition  of  A»  (1.4)  are  unknown,  probably,  do  not  exist. 
However,  some  sufficient  conditions  can  be  easily  established.  They 
are  given  in  Section  2,  In  Section  3  we  apply  them  to  matrices 
corresponding  to  symmetric  structures  of  finite  size.  Matrices 
associated  with  regular  structures  of  finite  size  are  quasi  block 
Toeplitz  matrices.  In  some  particular  cases  they  can  be  explicitly 
block  diagonalized.  Section  4.  We  propose  in  Section  5  a  special 
numerical  method  for  solving  system  (1.10)  with  quasi  block  Toeplitz 
matrix  k^  in  general  case,  Fast  Block  Elimination.  Finite  difference 
block  matrices  corresponding  to  the  biharmonic  operator  over  a 
rectangular  region  are  considered  in  Section  6.  The  asymptotic 
explicit  spectral  decomposition  of  such  matrices  is  constructed.  It 
leads  to  the  explicit  formula  for  natural  frequencies  of  a  clamped 
rectangular  plate  which  contains  the  Navier  formula  for  a  simply 
supported  plate,  as  a  particular  case. 

4.  All  proposed  decompositions  may  be  viewed  as  an  extension  of 
the  Poisson's  solver  [1,2].  Fast  Block  Elimination,  which  is 
considered  in  Section  5,  is  based  on  an  idea  of  simultaneous 
elimination  of  all  even  subvectors  X2>Xn,...  of  (1.10)  and  taken  from 
the  Fast  Fourier  Transform.  It  was  also  utilized  in  [1],  the  "cyclic 
odd-even  reduction  and  factorization  (CORF)".  However,  such  an  idea, 
being  applied  alone,  is  not  effective.  As  said  in  [1],  "from 
computational  viewpoint,  the  CORF  algorithm,  as  developed  here,  is 
virtually  useless."  We  combine  it  with  a  special  numerical  procedure 
for  computation  of  a  new  system  at  each  step.  This  makes  Fast  Block 
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Elimination  more  effective.  Its  efficiency  is  compared  with  the 
One-way  Dissection  Method  [3].  It  is  found  that  the  proposed  method  is 
more  effective  when  a  =  n/m  >  1.3  (see  Section  5). 
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2.   SUFFICIENT  CONDITIONS  FOR  EXPLICIT  BLOCK  DIAGONAL  DECOMPOSITION 
OF  MATRIX  A* 


1 .   THEOREM  2.1 : 

If  all  blocks  A,.  ^  of  A»  commute,  then  there  exists  the  explicit 
block  diagonal  decomposition  of  A^  (1.4). 


PROOF;   Suppose   all   A^.,   i,j=1 n  are  commutative  matrices: 

[a^j.A^  j  ]  =  Ojj,,  i.iiJ.Ji  =  1,...,n.  Then  all  A^j   have  the  same 
modal  matrix  V,  i.e., 


A.j  =  VN-jV"""  ,      i,j  =  1,...,n  (2.1) 


where 


Nij  =  r  Vij,  J';^^  ^   V  =  [y,X,.=^   ,   v-^  =  [%,]^,,=i  (2.2; 


Thus 


;.  A,  =  [VN^.V-Mjj^i  =  (I^xV)N,(I^xV-b  . 

Where  li^   is  a  full  block  matrix,  but  its  blocks  are  diagonal 

and  "x"   is  the  Kronecker  (tensor)  multiplication.  Therefore  one  can 
write 

A»  =  [(I^xV)P^]  [pIn^P^]  [P,(l^xv"b]  =  X^N^x;^       (2.^) 
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where  N»  is  block  diagonal  with  blocks  of  order  n 

N*  =  f  N,  J-^,  =  [  [v,jj?^.^,  J-^^  (2..) 

and  X#  is  a  full  block  matrix  with  rectangular  mxn  blocks 


THEOREM  2.2:  If  all  blocks  A^^^  of  A»  (1.5)  commute  and  therefore 
they  have  the  following  explicit  spectral  decomposition 


^T  =  "^ot""^  •      0'^=^ "^  (2.7) 


where 


"oT=  ^^Jo.  JJ.,  .    U=  [u,^]l^.,    .         0-1  =[uij]^.j.i  ,      (2.8) 

then  matrix  k^      (1.^)  has  the  following  explicit  block  diagonal 
decomposition: 


A^  =  U»M^U;1  ,  (2.9) 


with 


^*  =  i  ^J  J^l  =  ^  ^JctIo.t^I  Jj  =  1  •  (2.10) 


PROOF:  According  to  Theorem  2.1 


Therefore 


A*  =  P»A»pi  =  (P^Y^)M»(Y;'pi)  =  U,M»U;' 


with 


(2.11) 


A,^  =  Y^M^Y;^  ,  (2.12) 


where  M#  is  of  the  form  (2.10)  and 


Y»  =  (I„x{3)pT  ,      y;1  =  (I^xU"l)P^  (2.13) 


U,  =  P,(T  x{i)pl  -   U.I. 


u;i  =  p*(i^xO-^)pi  =   U~^Ij^ 


Note  that  in  (2.9)  matrices  XS^   and  M^  have  blocks  of  the  same  order  m 
as  matrix  A#. 

2.   THEOREM  2.3: 

If  both  matrices  A^^  ( 1  .  ^4 )  and  A#  (1.5)  contain  commutative  blocks 
then  A«  has  the  explicit  spectral  decomposition 


A»  =  w»  A,  w;i 


(2.1i4) 


Where  A»  is  the  spectral  matrix 


A.  =  [  [  A.   J-^,  J5^^  .  (2.15) 


n  m 

^JT  =   ^  ^stT^js^tj  =  ^    ^ja6\a^6T 
s,t=1  a, 6=1 


(2.16) 


and  W»  is  the  modal  matrix 


w.  .o.v.[[uijv„];^,.,];^j., 


(2.17) 


To  prove  this  theorem,  note  that  matrix  A»  now  possesses  two  block 
diagonal  decompositions  given  by  (2.3)  and  (2.9).  It  follows  from  here 
that 


M»  =  (U-^I^)(I^xV)N»(I^xv"'')(ijxl_^) 


n 

(0-^V)N,(UxV'b  =  [V[  I    N3^.Ui3U^.)V-M;,j=i 

s,t=1 


^^n^'^i    I         ^st^is"tj];,j  =  1^In'<^"') 
s,t=l 


(I^xV)A,(In'<V"'')  (2.18) 


Matrix  M^^  (2.10)  is  block  diagonal,  hence,  all  off-diagonal  blocks  of 
Ajt  are  zero  submatrices 


■10- 


n 

^    Ng^UigU^j  =  0„  .     t  *  s 
s,t=1 


Since  all  N  ^  (2.2)  are  diagonal,  so  are  the  diagonal  blocks  of  A» 


n 

im 


^j  =  ^    Ng^u.gU^j  =  [  Aj^  J^^^  .     j  =  1 n 

s,t=1 


I.e. , 


n 

^JT  =  ^    ^sti^is^tj  •     J=^ "  •   ^  =  ^ "^ 

s,t=1 


Thus, 


A^    =  (UxI^)M»(U''''xI^)  =  (U.<V)A^^(U"''xV"'') 


Finally   we  note  that  A^^   (1.5)  also  has  the  explicit  spectral 
decomposition 


A^^  =  W»A»W"''  ,  (2.19) 


where 


A.  =  [  A,  j;.,  =  [  r  Aj,  J"=,  \U    '  (2.20) 


m 

a,B  =  1 


W#  =  VxU  ,      W,t  =  V''''xu"''  (2.21  ) 
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3.  Consider  one  particular  case  of  matrix  A#  (1.'4)  which  is  the 
basis  for  next  sections: 

P 
A»  =   I   B^xC^  ,  (2.22) 

r=1 

where  B^  and  C^,  r=1 p,  are  nxn  and  mxm  matrices,  respectively. 

Since 

„  P  P 

A«  =  PiCI  B^xC^^P^  =  I     CpXBp  (2.23) 

r=1  r=1 

we  have 

P   (r)  -      P   (r) 

i,  j  =  l  ,  . . . ,n;    0, T=l , . . . ,m 


Thus,  one  can  formulate  the  following  statements: 


THEOREM  2,M: 


(i)  If 
B^  =  UM^iJ-l  ,    M^  =  [  Mj'^  JJ^i  .    U  =  [u.jlj^j^i  ,    r=1 p   (2.25) 


then  there  is  the  explicit  block  diagonal  decomposition  of  A»  (2.22) 

A»  =  U»M^jU;^  ,  (2.26) 

where  U>.   is  of  the  form  (2.11)  and 
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P  P   (r) 

M»  =   I  M^xC^  =  f   I   pj   C^  JJ^i  (2.27) 

r=l  r=1 


(ii)  If 


C^  =  VN^V-^  ,    N^  =  [  v/^  j;^^  .    V  =  [v^J-J^^^T  .     r  =  l p   (2.28) 


then 


A»  =  X^N^X;^  ,  (2.29) 


where  X^  is  of  the  form  (2.5)  and 


P  P   (r) 

N*  =   I  N^xB^  =  r   I  v^  B^  J!J^^  (2.30) 

r=1  r=1 


(iii)  If  conditions  (i)  and  (ii)  are  satisfied  simultaneously, 
then  A»  (2.22)  possesses  the  explicit  spectral  decomposition 


A^  =  W^A^,w;^  ,  (2.31) 


where  W^^  is  of  the  form  (2.17)  and 

P   (r)  (r) 


A*  =  r  r  I  Vj   v^   \U    JJ=,  (2.32) 

r=l 


REMARK;    Last   statement   is  well-known   [iJ] 


-13- 

3.   BLOCK  MATRICES  CORRESPONDING  TO  SYMMETRIC  STRUCTURES 
OF  FINITE  SIZE 


1.  We  define  symmetric  structure  (model)  of  finite  size,  S,  as  a 
physical  system  which  possesses  a  point-symmetry  group  G  of  order 
n  >  1.*  Clearly  each  such  structure  contains  n  identical  elementary 

regions  (or  cells)  S^ ,S2 S^.  We  call  S^   the  fundamental  region. 

Suppose  it  has  m  degrees  of  freedom,  then  matrix  A»  corresponding  to  S 
is  a  block  matrix  of  order  mn  and  can  be  presented  by  (l.M). 

2.  Constructing  matrix  A«  we  shall  obey  the  following  Symmetry 
Rule  which  states  that  we  are  free  to  choose  all  variables  and 
coordinate  system  for  the  fundamental  region  only.  The  variables  and 
coordinate  system  corresponding  to  region  S^  should  be  obtained  by 
applying  the  symmetry  transformation  g.  e  G,  j=1,...,n,  to  the 
fundamental  regions  S^ • 

3.  Let  gj^.g^  be  elements  of  the  group  G  and 

8w(i,j)  "  Sigj  .     i.J=1 n  (3.1) 

Introduce  the  following  n  permutation  matrices  of  order  n 

Q(gj)  =  [qik(gj)]?.k-l  =  t«k,w(i,j)];.k=1  '    J=1.---.n    (3.2) 

which  are  in  one-to-one  correspondence  with  the  group  table.  Denote  by 
T  (r=l,...,H)  the  r^^  irreducible  representation  of  group  G.  Its 
matrices  are  n  xn^,  unitary  matrices  (we  consider  finite  groups) 


*There  are  fourteen  point-symmetry  groups,  and,  hence,  fourteen 
different  types  of  symmetric  structures  of  finite  size. 


-■^H, 


r        IP 
TpCgj)  =  LTpag(gj)J^'^g^^  ,     r=1 H;  j  =  1,...,n 


(3.3) 


For  point-symmetry  groups  1  <  rip  £  5.  Elements  TpctB^^i^  °^  ^'^^ 
irreducible  representations  satisfy  some  orthogonality  relations,  in 
particular, 


H  n^  "r 

^  ~  ^    "^raB^^i^VaB^^k^ 
r=1    a,B=1 


=  6,- 


ik  ' 


i,k=1  , . . .  ,n 


(3.M 


LEMMA  3- 1 :  There  exists  the  explicit  block  diagonal  decomposition  of 
matrices  Q(gj)  (3-2) 


Q(gj)  =  UT(gj)u"  ,    j=1 n 


(3.5) 


where  T(g^)  are  block  diagonal  unitary  matrices  having  the  same 

H    — T 

configuration  for  all  g^  £  G,    and  U   is     a     uninormal     matrix,      U     -  U    , 
namely, 


T^^^)  = 


Aj    TVVAE© 


X 


H%^ 


"^^^v^ 


J'-J 


rv^    TiMeS 


(3.6) 
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or  in  short 


T(gj)  =  [  r^,v(gj)JT=1  Jr=1 

=  f  F  [SYVa6(gj)]I'B=iJYliJ"=i   '   J  =  l "   (3.7: 


and 


■"ni/2rr.   r.  ^l"r  i"r  in.H 


"=Uis];.s=1  =  U-f)^^'[[VTa(Bi)Ul  V=^    ]?:?=1        (3.8) 


Here  column  subscript,  s,  is  associated  with  three  subscripts  r,  a  and 
Y.  Since  ^^^.-.n^  =  n,  subscript  s  runs  from  1  to  n  when  a  and  y  run 
from  1  to  n  and  r  runs  from  1  to  H. 


""  =  l=sk!s",k.,  -  [(T)"'[["r.6(8.>l6M  ]":A".,  '3-9' 


PROOF;  Equation  (3,5)  is  a  matrix  form  of  the  following  orthogonality 
relation 

H  n     "r> 

I     -^         I  VYa^Si^VaB^Sj^VYB^gk)  =  «k,w(i.j)  •      (3-10) 

r=1    a,6,Y=1 

i.j.k=l, . . .,n 
which  generalizes  (3.'^). 


i».   THEOREM  3.1: 

Block  matrix  A«  (1.4),  associated  with  symmetric  structure  S  and 
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constructed  in  accordance  with  the  symmetry  rule,  can  be  expressed  by 
the  following  Structural  Formula: 

n 

A»  =   I  Q(gj)xA^j  ,  (3.11) 

J  =  1 


where  A^  ^  are  blocks  of  the  first  block  row  of  A^ 


PROOF:   Introduce   the  identical  symmetric  structure  S'  —  S  with 

elementary  regions  S,  =  S ■ ,  S,  =  S..,    .  Suppose  that  we  coincide  S 

I    J   <--         J  I  »  •  •  • 

» 
and  S  by  a  certain  symmetry  transformation  g  £  G  so  that  Si  be 

coincided  with  S.  .  Then  Si  "iH  coincide  with  S^,.     .x,  where  w(i,j)  is 

defined  by   (3.1).   Suppose  matrix  A^^  =  [a  ■]"  ^_,  corresponds  to  S 

PJ   P » J  - 1 

while  A»   =   [Apj]pj  =  i    to  s'.     Clearly 


Vj    =  ^v(i.p),w(i.j)    •  i'J'P=l " 


On  the  other  hand,  S  and  S  are  identical,  thus,  their  matrices  while 
satisfying  the  symmetry  rule,  have  to  be  equal.  Hence  A^  =  A  . . 
Therefore 


^pj  -  ^w(i,p),w(i,j)  '     i,j,p-1,...,n 


Let  p  be  the  equal  to  one.   Then,   since   w( i, i)  =  i,    (3.1). 
^j  =  ^i,w(i,j)  o-^ 


^1j  =  ^  Ai)l^il,w(i,j)  •    j  =  1....,n 
il=1 
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Multiplying  both  sides  by  6,^  ^/^    .^   and  summing  up  with  respect  to  J 
from  1  to  n,  we  obtain 

n 

^  '^ij^k.wCi.j)  =  ^ik  •    ^'^='^ " 

J  =  1 

or,  by  virtue  of  (3-2) 

n 
^ik  "  ^  ^ik^Sj^'^lj  '    i,k=l,,,.,n  (3-12) 

j  =  1 

the  Structural  Formula  (3-11)  follows  from  here. 


5.   THEOREM  3-2: 

Block  matrix  A«  (3.11)  has  the  following  explicit  block  diagonal 
decomposition: 

A»  =  U^A^^U^  .  (3.13) 

where  U^   is  uninormal 

U*  =  UxIn,=  [Ui3lj?.s=l  ^3.1^) 

and  ^y^   is  Hermitian,  block  diagonal 


n 


'r   |H 


A*  =  I  T(gj)xA^.  =  [  [  6,/,  J,:^  j;!.i  (3.15) 

j  =  1 


with  blocks  A^  of  order  mn^,  1  <.  n^,  <  5 
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n  n  n 

J  =  1  j  =  1 


PROOF:   Equations  (3 . 1 3)-(3 . 16)   follow  from  (3.11)  if  we  substitute 

there  (3.5)-(3.9). 

a 

NOTE:  Symmetric  structure  S  has  at  least  n^-fold  natural  frequencies 
and  critical  loads,  1  <.  n^  £  5,  because,  as  it  follows  from  (3-15) , 
matrix  A^^  has  at  least  n^-fold   eigenvalues. 


THEOREM  3.3:  Block  matrix  A»  can  be  presented  in  the  form 


n 


a;^  =  I     Q(gj)xA^j  .  (3.17) 

J=1 


where 


H  n   "r    _   _ 
r=1    a, 6=1 


and  A„  o  are  blocks  of 
ra3 


I'V 


K'  -  [^^6=1  ^'-''^ 


PROOF:  We  find  from  (3.11)  that 


n 

a;^  =  I  Q(gj)xA^j  , 
j  =  1 
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where  A^ .  are  blocks  of  the  first  block  row  of  A».  On  the  other  hand 


and  therefore 

"  _   n 

a;^  =  I  T(gj).Aij  =  [  r  SyV^  Jy=i  J?=i 

Suppose  blocks  A~\  r=1,...,H  are  computed.  Then  according  to  (3.16) 

_  n 

VaB   =     I      VaB^8j)Aij    ,  a, 6=1 n^;      r  =  1,...,H 

j  =  1 

Multiplying  both  sides  by  (np/n)Tj,^g(g|^)  and  summing  up  with  respect  to 
a  and  6  from  1  to  n  and  r  from  1  to  H,  we  obtain,  taking  into  account 
the  orthogonality  relation  iS-^), 


r=1    a, 6=1 

nHn"r  _  n_ 

j  =  1     r=1    a, 6=1  j  =  1 


Thus,  to  invert  mnxmn  matrix  A»  (3-11)  we  need  to  invert  only  the 
matrices  A^  (r=l,...,H)  of  order  mn^,  1  <  rij.  <  5.  Note  that  if  A»  is 
symmetric,  some  blocks  A^ .  are  mutually  transposed,  we  have  to  compute 
only  truly  different  A^.. 
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6.  THEOREM  3.4;  Linear  system 

A»x»  =  b,  (3-20) 

of  order  tnn  is  split  into  H  decoupled  subsystems  of  order  mn^, 

1  £  n^  <^  5,  containing  nj,  unknown  subvectors  Ypy*  ^~^ "p'  ^^^^ 

Ar[ypi  .....ypn  ]  =  [^r] '^rn  ^  •  r=1,.,.,H          (3.21) 


where 


rY  =  [^V^^    [    I     bkT^Yg(gk)]6=i    Y=1,....n^;  r  =  1 H   (3.22) 

k=l 


c   -  '  '^^1/2 


Once  (3.21)   are  solved,   the  original  unknowns  x»  =  [xj^]"_i   are 
determined  by 


H   n 

r)^/2 


Xi  =  I      (-J-)     I  yrYaVYa^Si)  ,      i  =  l....,n    (3.23) 
r=1         Y,a=1 


where  YpYa  ^^^  subvectors  of  y^^^  ^  l-yrYa-la=1*  ''""  ^='*^^»  substituting 
(3.13)  into  (3.20)  and,  taking  into  account  (3.8)-(3.9),  we  obtain 
(3.21)-(3.23). 

7.  Finally,  it  is  necessary  to  note  that  the  Structural  Formula 
(3.11)  can  significantly  simplify  the  use  of  tetrahedrons  and  cubes  as 
finite  elements  because  the  corresponding  matrix  A^^  is  symmetric  and 
therefore  it  is  sufficient  to  compute  and  store  only  one-half  of  its 
first  block  row. 
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4.   BLOCK  MATRICES  CORRESPONDING  TO  REGULAR  STRUCTURE  OF  FINITE  SIZE 

1.  We  define  the  regular  structure  (model)  of  finite  size  as  a 
structure  (model)  which  is  formed  by  a  finite  number  of  identical  (or 
periodically  repeated)  substructures  (submodels)  of  finite  size. 
Matrices  associated  with  regular  structures  are  quasi  block  Toeplitz 
matrices.  Their  blocks  at  the  left  top  and  right  bottom  corners 
express  the  boundary  conditions  and  therefore  they  usually  distinguish 
from  others  in  the  same  block  codiagonal. 

2.  Introduce  the  following  auxiliary  quasi  Toeplitz  matrix  of 
order  n 


\   = 


11  ^ 

\     0  \ 
\     0  1 


\   O  1 


r\ 


U\.\^\^  ^ 


(A.l) 


LEMMA  4.1:  There  is  the  explicit  spectral  decomposition  of  B^ 


B^  =  UM^U^  , 


(4.2) 


where 


(1) 


■  r  "J   ij 


(1) 

y^    =2  cos  6^ 


(4.3) 
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"=[-ij]lj  =  1  ^'-'^ 


u^.  =  c.(sin  ie.  +  f  sin  (i-1)e.] 

=  c.[g  sin  (n-i)e.  +  sin  (n+1-i)ej         {^.5) 

J  J  J 


Arguments  6.,  j=l,...,n  are  zeros  of  the  characteristic  polynomial 


A(e)  =  sin  (n+1)e  +  (f+g)  sin  ne  +  fg  sin  (n-1 )e      (4.6) 


There  are  n  distinct  zeros  6^  on  (0, it].  We  omit  the  proof  because  it 
is  the  same  as  for  Lemma  6.1. 

Introduce  a  Chebyshev  family  of  matrices  B 


Bo  =  In 


Bl 

B2  =  B^  -  2Bq  (4.7) 

B3  =  b3  -  3B^ 

B^  =  of  -  4b2  .  2B^ 


They  are  quasi  Toeplitz  matrices  and  are  given  explicitly  by  (4.8)  thru 
(4.11)  for  n=9. 
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1.   = 


P^-^^^ 

-"s 

K 
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^      '1            1 
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K- 


Rl^-^'V-(^-^'^ 

-h 

^ 

-e^-^^^ 

-^ 

■  -  —  - 
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Clearly, 


B^  =  UM^U' 


(M.12) 


with  the  eigenvalues 


(r) 
y.   =2  COS  re.  , 


j  =  1 n 


(M.13) 


Introduce  also  a  family  of  non-symmetric  matrices  B^  of  the  same 
order  n 


^o  "  ^n  •    ^1  • 


B2  =  b2  -  2aBQ  , 


i^.^^) 


Their  explicit  forms  are  given  by  {^  .^5)-{^  .^6)   for  n=9 


A 

i. 


"i  '  '^ "' 

"   1  Po. 


I        J 


(4.15) 


-\ 

\ 

\ 

""] 

oC^ 

f. 

V 

\ 

d^"- 

-'^i      ' 

"^'^~~' 



oci 

-? 

vt.^ 

.^CJ 

-^'^-ri 

(4.16) 
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LEMMA  4.2:  There  exists  the  explicit  spectral  decomposition  of  B_: 


B^   =   UM^U    '     , 


(4.17) 


where 


y.        =   2a  cos  re^    , 


i/2 

^ij   =  «         "ij    ' 


"ij   =  " 


■J/2 


CJ.IS) 


(4.19) 


'ji    • 


This  lemma  immediately  follows  from  Lemma  ^4.1  if  we  note  that  there  is 
a  diagonal  matrix  S 


s  =  [  aJ/^  j; 


J/2  in 


(M.20) 


such  that  S  B  S  is  symmetric  quasi  Toeplitz.  For  instance. 


'p'^,   ^  = 


-1- 

O 

f^ 

f^ 

/^ 

o 

v]^ 

^ 

Then 
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u  =  su  , 


(4.21) 


where  U  is  of  the  form  (4.4). 

3.  Consider  the  following  two  quasi  block  Toeplitz  matrices  of 
order  mn 


r=0 


p  <  n-1 


(4.22) 


and 


r=0 


p  <  n-1 


(4.23) 


where  A^  (r=0, 1  , . . . ,p)  are  arbitrary  mxm  matrices.  They  are  given 
explicitly  by  (4.24)  for  n=7,  p=3  and  by  (4.25)  for  n=7,  p=2, 
respectively. 


k.- 


K-SK 

A. 
A, 

A. 

A, 

K'\k 

A. 

A, 

A, 

A. 

K 

k. 

A, 

A, 

A. 

A, 

A. 

k 

k. 

Az 

A. 

A, 

A, 

Ai-^^i 

A. 

A. 

A, 

Ao-?  Aj 

A. 

K-'ik 
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k,- 


^^-hK 

K 

^(>^^^^ 

k. 

k. 

k. 

oC^A, 

«L^, 

k. 

A^ 

A. 

^"-K 

oi  ^v 

K. 

A, 

^. 

Ji'k^ 

^  A\ 

A. 

A. 

^. 

ocH, 

o(  A^ 

A. 

Ar^^i 

o^^A, 

4Ar^^^ 

\S\,\%\i\ 


THEOREM  4.1:  Matrices  A»  (4.22)  and  A«  (4.23)  have  the  explicit 
block  diagonal  decompositions 


U»T»uI 


(4.26) 


and 


A«   =  U»T^tU, 


(4.27) 


Here 


U»   =  UxT      ,  U»   =  U»   =  Uxl.   =  SUxI 


(4.28) 


T»   =      I     M^xA^    =    [  Tj    JJ^T 
r=0 


(4.29) 
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r=0 


P 

Tj  =  Aq  +  2   I  A^  cos  rGj  ,  (^.31) 

r=1 


P 

Tj  =  Aq  +  2   I  A^a'^^2  cos  j,g^  ^  (M.32) 

r=1 


THEOREM  4 .2:   If  matrices  Aj, ,  r=0,l,...,p  commute,  i.e.,  they  can 
be  presented  in  the  form 

Ap  =  VN^V-1  ,  N^  =  r  v/  J'J^T  ,     r  =  0,1 p      (4.33) 

then  Aj(   (4.22)  and  h^      (4.23)  possess   the   explicit   spectral 
decompositions 


A»   =  W,A»w|  (4.34) 


and 


A,  =  W»AjtW»^    ,  (4.35) 


where 


w,j  =  uxv  ,         w^f  =  Uxv  =  suxv  ,  (4.36) 

A*    =    [    r   ^u   Jx  =  l    Ji=l     .  A,    =    r    [   X.      J-^,    J5^^     .  (4.37) 


(o)  P        (r)  .  (o)  P        (r)      ^/2  ,,      „^ 

Xj^   =   v^        +21      v^        cos  r6j    ,        Xj^   =   v^        +   2      1      v^        a"^ '  "^  cos   re^      (4.38) 

r=1  r=1 
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Clearly,  both  theorems  are  particular  oases  of  Theorem  2.M.  Symmetric 
quasi  block  Toeplitz  matrices  occur  in  applications  usually  more  often. 
Note  that  matrix  A#  (4.22)  becomes  symmetric  if  all  submatrices  A  are 
symmetric.  Matrix  A^^  {^.23)  will  be  symmetric  if  (i)  all  even  blocks 
Aq,A2,...  be  symmetric,  (ii)  all  odd  blocks  A^.A^,...  be 
skew-symmetric,  and  (iii)  a  =  -1  ,  f  =  g  =  o. 

^.     Consider  Hermitian  quasi  block  Toeplitz  matrices.   With  no 
loss  in  generality  consider  tri-block  diagonal  matrix 


K- 


A',    A.   A 


a:  k.  k 


A  O       ^ 


a:  a" 


mn 


(4.39) 


THEOREM  4.3:  If  all  blocks  of  A»  (4.39)  commute,  i.e.,  if  they  can 
be  expressed  as 


Aq  =  UMU^ 


UNU' 


Aq  =  um'u" 


A*,^  =  UNU" 


II       »  '  H 
Aq   =  UM  U" 


(4.40) 


where 


U  =  l%,r,,,-.,     ,  M  =  [  u,  \%     .     N  =  [  v^  J^^^  . 


(4.41) 
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then  there  exists  the  following  explicit  block  diagonal  decomposition 
of  A»  (i).39): 


where 


A»  -  X^TjfX,   , 


(M.iJ2) 


T 


(M.M3) 


(14.44) 


T.  = 


(4.45) 


r\ 


This  theorem  is  a  particular  case  of  Theorem  1.1.   Clearly,  if 


|fj,|gj    <   1    .  1=1 t 


(4. 46) 


where 
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N"^  ^"^ 


f,  =-ri-f.     g^=_L^.     t=1 m         (M.M7) 


|vj         ^     |v^ 


then  A,   (M.39)  has  the  explicit  spectral  decomposition  (Theorem  1.3). 
It  will  be  also  true  if  |f|,|g|  <   1,  where 

a;  =  A^  -  f(u|N|u")  .     a;'  =  A„  -  g(U|N|u")  .     |n|  ^  [  |vj  J?^,  i^.^8) 
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5.   FAST  BLOCK  ELIMINATION  FOR  SOLVING  QUASI  BLOCK  TOEPLITZ 
SYSTEMS  OF  LINEAR  EQUATIONS 


1.  Without  loss  of  generality  consider  the  following  symmetric 
tri-block  diagonal  system  of  order  mn  with  blocks  of  order  m 


A'    'h 


2>^  A  ^ 


h. 

^. 

Xz 

K 

X. 

= 

^. 

^n-\ 

W-\ 

i. 

V.n 

(5.1) 


»    »  ? 


Matrix  A»  contains  only  four  different  nonzero  blocks  A,  A  ,  A  ,  and 
B.  All  numerical  methods  [3]  which  can  be  applied  to  (5.1)  ignore  this 
property,  moreover,  they  totally  destroy  it  during  elimination  (block 
elimination).  However,  we  can  preserve  and  exploit  this  essential 
property  if  at  each  step  of  the  forward  pass  we  shall  simultaneously 

eliminate  all  even  subvectors  X2,Xj^,x^ Then  a  new  system, 

obtained  at  each  step  of  such  block  elimination,  will  have  the  same 
block  periodic  configuration  (5.1),  but  its  order  will  be  approximately 
half  the  previous  one: 


n^^^  =  E(n^/2)  =  e[  (n^-l  )/2'<^^  ]  +  1  , 


"o  =  "  ' 


k=0,1,2,...   (5.2) 


E(x)  is  the  greatest  integer  of  x. 
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2.  This  idea  is  very  attractive,  however,  it  becomes  efficient 
only  if  we  shall  accompany  it  with  a  special  numerical  procedure  in 
computation  of  desirable  blocks  Aj^+^  ,Aj^^^  ,  A^^^^  ,  B,^^.^  and  subvectors 
h^*\  k=l,2,...  at  each  step.  Let  us  compute  them  from  the  following 
equations 


Vl  =  A,  *  U  .  W 


Vl  =  ^ 


Vi  =  ^  *  " 


r\    *  w 


'k+1 


for  odd  n, 


-A|^  +  U  +  W    for  even  n^^ 


(5.3) 


b'<^^  =  b!^.  . 


^k.1 
"k+1 


2j    ^2j    ^2j-2  ' 


\  '  %-i 


\-^ '  \  •  '\-^ 


j  =  1,2,...,n|^^^  -  1 


for  odd  n. 


for  even  n, 


(5.4) 


where  blocks  U,  V,  W,  U  and  subvectors  Zp^  ,  y2i-2'"*  ^^^  found  by 
standard  Gaussian  elimination  applied  to  the  following  augmented 
matrices 
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A. 

il 

Ik 

X' 

X 

o 

o 

...0... 

...0... 

^l 

o 

,^Ki  S.    '     T 


C. 


V'  I  Atr  \-i 


(5.5) 


and 


XI:  i: 


\  iO 


-- 

1 « J 

— p-* 

~1-,^K 

0 

u 

^ 

^"«_ 

(5.6) 


The  last  transformation  has  to  be  done  only  if  n^^  is  even.  Certainly, 
only  upper  triangular  parts  of  U,  W  and  U  are  computed.  After  p  steps 
of  such  block  elimination,  where 


p  =  E(logp  (n-1))  , 


(5.7) 
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we  obtain   the  system  of  order  mrip  =  2m 


P 


P 
PJ 


xP 


L    ^J 


^'2 


(5.8) 


which  does  not  contain  repeated  blocks  and  has  to  be  solved  by  any 
standard  method. 

3.  Back  substitution  also  requires  p  steps.  At  the  k  step  of 
the  backward  pass  we  compute  subvectors  xR  ,x^  ,...  which  have  been 
eliminated  at  the  p-k  step  of  the  forward  pass.  Upper  triangular 
matrices  Rjj-k  (5.5)  and  R  _,  (5.6)  are  used  in  this  case: 


Y-k 


[Xp   fXjj   ,...,Xj.    1    ~    1^2      •''1)   ••••iC^    J 


'p-k 


p-k 


'■^p-k'^p-k- 


,p-k+1 


p-k+l 
X2 


D-k+1 


p-kH 

p-k+l 
^p-k+1 


(5.9) 


where 


'p-k 


p-k 


,   if  n  1^  is  odd 


n   _-,  ,   otherwise 
p-k  ^ 


(5.10) 


4.  Clearly,  each  step  of  this  procedure  requires  more  computation 
than  standard  one,  but  a  number  of  steps  is  very  small.  Therefore  the 
efficiency  of  Fast  Block  Elimination  will  increase  with  a  =  n/m 
increase. 
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Compare  this  method  with  the  One-Way  Dissection  Method  (1WD)  [3].* 
Let  (N^,M^)  be  a  numer  of  operations  and  storage  requirement  of  the 
proposed  method,  respectively.  The  same  quantities  for  the  1WD  are 
(Np.Mp)  and  (Np  .Mp  ),  where  N2  and  Mp  are  optimal  number  of 
operations  and  optimal  storage  requirement,  respectively.  Then  the 
efficiency  of  Fast  Block  Elimination  can  be  characterized  by  the 
following  four  ratios  as  functions  in  m  and  n 


and 


M**  M** 

»»    ''2  **     2 


Corresponding  curves  R(m,n)  are  depicted  in  Figures  1  and  2.   From  them 
it  is  concluded  that  the  proposed  method  is  preferable  (R  >  1 )  when 


*The   Nested   Dissection  Method   (ND)   is  more  sophisticated. 
However,  it  has  approximately  the  same  efficiency  as  the  1WD  [3]. 
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a  =  n/m  >  1 ,3*  (5.11) 

5.  The  proposed  method  is  also  applicable  to  the  eigenvalue 
problem.  Suppose  we  use  the  Sturm's  method  whose  procedure  requires 
factorization,  by  Gaussian  elimination,  of  the  matrix  B»(X)  =  AI-A»  for 
some  particular  value  of  X  and  to  compute  the  negative  index  of  inertia 
0  (i.e.,  a  number  of  negative  Gaussian  pivots  of  B»(A)  transformed  to 
a  triangular  form).  Then  applying  the  Fast  Block  Elimination  one  can 
take  advantage  of  repeated  blocks  in  computation  of  o  . 


*More  precisely:  if  n  <  40,  a  >  1.5,  if  n  =  50-60,   a  >  1.1-1.2,   if 
n  >  70,  a  >  1 . 
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EXPLICIT  ASYMPTOTIC  SPECTRAL  DECOMPOSITION  OF  FINITE  DIFFERENCE 
BLOCK  MATRIX  CORRESPONDING  TO  BIHARMONIC  OPERATOR  OF  A  REGULAR 
RECTANGULAR  MESH 


1 .  The  1 3-point  molecule  for  a  biharmonic  operator  on  a  rectangle 


A.  ^y  = 


't   "X 


C.xV 


1    r 

^        S 

-t 

\ 

S 

? 

i 

A 

t 

^ 

t 

y 

(6.1) 


with  boundary  conditions 


edge 


=  0 


and 


9n  ledge 


or 


92u 


9n 


2  I  edge 


leads  to  the  following  five  block  diagonal  quasi  Toeplitz  matrix  A,  of 
order  mn  with  blocks  of  order  m: 
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Im 

2> 

Im 

:5> 

k 

i. 

T. 

A.= 

•    •    . 

•         •> 

•          •          • 

»                          • 

• 

X^^ 

:b 

K 

2> 

In. 

Im 

:fe. 

A^Io. 

(6.2) 


where 


A  = 


>ir 


^      ^ 


^        f      ^       ^ 


r      c^      o 

o^     r 

r 

^  ?^' 


2p- 


^     -4: 

^     S   4: 


4:      S     4 


(6.3) 
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p  =  6  +  8y^  +  ey^ 


q  =  -i4T^(1+Y^) 


r  =  Y 


s  = 


■4(1+Y^) 


2Y' 


(6.4) 


Y  =  ^ 

Ay 


Ax 


X 

rviT 


Ay 


tn+1 


(6.5) 


Matrix  A»  (5.2)  and  all  its  blocks  are  symmetric.  The  boundary  blocks 
of  A»  (6.2)  and  boundary  elements  of  A   (6.3)  are  A  ±  I^j,  and  p  ±  r, 

respectively.   Positive  sign  corresponds  to  8u/3n  =  0  (clamped  size), 

2     2 
while  negative  sign  to  9  u/9n  =  0  (simply  supported  size).   If 

boundary  blocks  of  A»  are  A-I^^  or  if  boundary  elements  of  A  (6.3)  are 

p-r,  matrix  A*  can  be  explicitly  block  diagonalized.  If  both  such 

conditions  are  satisfied  (a  simply  supported  rectangular  plate),  matrix 

Pi^   has  the  explicit  spectral  decomposition,  which  can  be  viewed  as  a 

discrete  Navier  solution.  We  shall  consider  a  case  with  A+I^^j  and  p+r, 

i.e.,  a  clamped  rectangular  plate. 

2.  In  accordance  with  Lemma  4.1  matrix  B  (6.3)  has  the  explicit 

spectral  decomposition 


B  =  VMV  , 


^-\     ^.  lU    ' 


s  +  2t  cos 


Tir 
m+1 


(6.6) 
(6.7) 
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V-1  =  V  =  [v„  1";  ,  ,  .    v„  =  /_L  sin  -^        (6.8) 


LEMMA   6.1:   There   exists   the  following  explicit  spectral 
decomposition  of  matrix  A  (6.3) 

A  =  UAU'^  ,  (6.9) 

A  =  [  A^  J'^^^  ,     A^  =  p  +  2q  cos  6^  +  2r   cos  2rQ^    ,     (6.10) 
U=  [u„J^.,  =  i  .    U"^  =UT  .  (6.11) 


cosh  (m+1  )4<^  -  cos  (m+1)e^   sinh  o^^       sin  06^ 

u_^  =  c^[cosh  04)^  -  cos   a6^  -  ( -  1] 

°^  ^  ^  ^         sinh   (m+1)ii)         sin    (m+1  )6        sinh  \\)  sin   9 


sinh  ^^  sin  9 


(6.12) 


with  \\)^   determined  from 


cosh  i|;  =  -  -3-  -  cos  9^  ,     T=1 m         (6.13) 

T     2r        ^ 


Arguments  6^  are  zeros  of  the  following  characteristic  polynomial 


•    2  -2 

A^(9)   =   2[l    -  COS   (m+1)9   COS    (m+Dti;  -  — — ^  "  ^^""^     '^   sin    (m+1)9   sinh   (m+1)i|;] 
'"  sin   9  sinh  \li 

=   [sinh  \\)  tan  — —  9  -  sin  9   tanh  ti^J    [sin   9   tan  9   +  snh  i|;  tanh  \li 

(6.1i4) 
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where  \\)   is  defined  by 


cosh  \\)  =  -  -L  -   cos  e 
2r 


(6.15) 


PROOF;  Present  equation  (A-AIj^)u  =  0  as  the  following  finite 
difference  system 


ru 


0-2  ■"  ^Vl  *  (P-^)%  ""  ^%+l  ""  "^Vl  =  ° 


(6.16) 


u_^  -  u^  =  0 


u^  =  0 


m+l 


"m+2  -  "m  =  0 


(6.17) 


Present  u   as  n^  and  denote  y  =  n  +  n  •  Then  we  obtain  from  (5.16) 

2 

'"y     +qy+p-2r-A  =  0.        Present     also     X      in     the     following     form 

A  =  p  +  2q  cos  e  +  2r  cos  26,  then 


y^  =  2  cos  9  ,    y2  = 


-  -  2  cos 
r 


=  2  cosh  4)  , 


^     is     real,      since  -q/r   =   4(1+Y  ^)   >_  4.      Thus,    n^    2   "   ^*^^  >    ^34   =   ^'^ 
and 


u     =   Ci    cos   c9   +   Co   sin   oG  +   c-,   cosh   oil)   +   c .   sinh  ail; 
0  1  2  34 


Substituting  into  (6,17),  we  obtain  (6.12)  and  (6.14) 


It  is  important  to  note  that  there  are  m  distinct  zeros  6^  on 
(O.tt] 


JL  <  e  <  iTllll  .      1=1, ...,m  (6.18) 

m+1    '^  m+1 


Finally,  we  present  6^  in  the  form 


(T  +  e^)TT 


m+1 


0   <  e^   <.5    ,  T=1 ,...,m  (6.19) 


3.      LEMMA   6.2: 


k^    =   I^xA    +   W*N»W|    ,  (6.20) 


where 


V^   =  Vj^Ujf  ,  (6.21  ) 

V»   =  I^xV  ,  V    (6.8)    ,  (6.22) 

">  =  Kj]?.j  =  i=[    \^  \Un,,  =  ^    '  ^6-23) 

N*    =    [   Nj  JJ.,    -    \     \   v]'^    J-^,    \-__,  .  (6.2i4) 

(t)  (t)  ^^  (O  (t) 

V.        =   2s  cos   Oi        +    4t   cos  cos  6.        +2   cos   20.  (6.25) 

J  J  m+1  J  J 
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(t)  (t)    r  (t)  (^) 

u..      =  C::        [cosh   ii|;^        -  COS   ie^ 


(x)  (t)  (t)  (t) 

cosh   (n+1)il..        -  cos   (n+l)e.  sinh  li|j-  sin   ie. 

J  J [ ^  -  1_)]   (6.26) 


(t)  (t)  '      (t:         (t) 

sinh  (n+1)ii).     sin  (n+1)e.     sinh  i|j.      sin  6. 

To  "      To 

sinh  \\i.  sin  6^ 


(t) 
where  4;^        are  found  from 


("^^  1  TIT  ^"f)  ,r         X 

cosh   il;-        =  -  -  s  -   t   cos  -  cos   9^  (6.27) 

J  2  m+1  J 


(t) 
Arguments  6^   ,  j=1,...,n,  are  zeros  of  the  polynomial 


(t)  r  Ct")  1 

A^      (6)   =   2[l    -  cos    (n+1  )e   cosh   (n+1  )i1j^  ' -^  ■• 


sin^  9  -  sinh^^f^^^   ^^^    (^.^  )g   3^^^,  (n^D^.^^^ 
sin   9   sinh  ^i^^'^ 


=  [sinh  \li^'^^   tan  I^  6  -  sin  9  tanh  ^  ip^"^^] 


[sin  9  tan  I^  9  +  sinh  i^^^^  tanh  I^  4;^^^]  ,   (6.28) 


where  tj;^  -^  is  defined  by 


cosh  il;^'^^  =  -  1  s  -  t  cos  -i^  -  cos  9  (6.29) 

2  m+1 
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To  prove  this  lemma,   present  block  matrix  A#   (6.2)  in  the 

following  form 


A,  =  I„-A  *    (I  xV)C*(I  xV)  , 


n   '  *  n 


(6.30) 


where  C»  is  a  block  matrix  with  diagonal  blocks 


c- 


K\  O  VA  T^ 
I^  M  O  NA  1, 


1^  M  O  VA  T^ 
I^  VA  O  VA 


(6.31) 


Then  matrix  P«C,P,t  is  block  diagonal  with  nxn  blocks  C^ 


P#CjtP^  -  I  C^  \x  =  -\    » 


(6.32) 


^x  = 


>  t=>\,,..,m 


(6.33) 
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Blocks  C  have  the  same  configuration  as  A  (6.3),  therefore  Lemma  6.1 
can  be  applied: 


C^  =  U^N^U^  ,      T=l m  (6.3i^) 

U,  =  [Uij  ]l^__,    ,  N^  =  [  Vj   J5^i         (6.35) 

(t)      (t) 
Elements  v.   and  u^.   are  determined  by  (6.25)  thru  (6.29).   Clearly, 

the  explicit  spectral  decomposition  of  C^^  (6.31)  is 


C,  =  P»  [  c^  Jpi 

=  [P.  [0,  J'^.^pI]  [P»f,  N^  J'^^iFI)  [P*f  OT  J'J^^fI) 

=  UjtN»U»  ,  (6.36) 

where  U^^  and  N»  are  given  by  (6.23)  and  (6. 24),  respectively.  Finaly, 
substituting  (6.36)  into  (6.30)  we  obtain  (6.20)-(6.22) .  This 
completes  the  proof. 


(t) 
REMARK:  Each  characteristic  polynomial  A^   (6),   T-1,...,m  has  n 

distinct  zeros  6.   ,  j  =  1 n  on  (O.ttJ 

IL  <  e^'^  <  ii^:^  .    J=1 n  (6.37) 

n+1    J      n+1 
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We  present  them  in  the  form 


(t)         [j+e-i      U  (t) 

'i        =  —, .  0   <   e.        <.5    ,  j  =  1 n;    t  =  1 m      (6.38) 

J  n+ 1  ^ 


It  is  important  to  note  that  there  exists  the  limit 


(t) 
Ej  =  lim  ej    ,  (6.39) 


which  is  the  same  for  all  T=1,...,m. 

^.   LEMMA  6.3t   Block  matrix  A»   (6.2)  can  be  presented  in  the 
following  form 


A»  =  X^H^xl   +  0(T^  ,  (6.ilO) 


where  H#  is  block  diagonal 


H»  =  [  Hj  JJ^^  ,      "j  =  '^  *  ™jV  (6-^1) 


and 


X^^  =  V»U^V»  (6.i^2) 
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PROOF ;  Introduce  n  matrices  X^j  similar  to  U^ 


r  (ij^m 


(6.M3) 


Their  elements 


(ij) 


2        V       ^°'^    -.._    oaiT    „,„   a  TIT 


OT  tn+1 


)      u.  .        Sin  — -  sin  — - 

m+1  iJ  m+1  m+1 


a=1 


(6.44) 


possess  periodic  properties  which  leads  to 


\}'\   = 


2r 


o 

^x    •  •  • 

o 

* 

^ 

D.. 

O 

0 

(6.45) 


i;^=  A,...,  n 


(ij) 


Since  r  =  Y""  and  x^^^       =  0(1), 


AX.j    =  X^jA   +   O(T^)    ,  i,j=1,...,n 


and 


AVU^j    =   VU^j(VAV)    +    O(Y^)    ,  i,j=1 n 


CI 


early,    VU-  •    =  W^^.    -  blocks  of  W»    (5.21).      Thus 
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AW.j  =  W.j(VAV)  +  0(y'')  ,      i,j  =  1,....n  (6. 46) 


or 


(I„xA)W,t  =  W,((I„xVAV)  +  O(Y^) 


n   '  *    *  n 


and 


I^xA  =  W»(Ij^xVAV)w|  +  O(T^)  (6.i47) 

Introduce  this  equation  into  (6.20) 

Obviously  we  can  write  it  in  the  form  of  (6.40). 

D 


(t) 
5.   Elements  v.    (6.25)  of  blocks  N.  (6.24)  can  be  written  as 


(t)  Ij+ei  U  2ij+e.   Jtt 

v^   =  2y^  cos  ^ +  2  cos ,    j  =  1,...,n;  T  =  1,...,m 

J       ^  n+1  n+1 


Therefore  according  to  (6.39) 


(t)  (j+eJi:  2(j+e.)iT 

lim  V-        =   2u      cos —   +   2   cos —  =   2y^   cos   9.    +   2  cos   26- 

J  T                n+1  n+1                   ^              J                         J 
n->-oo 

j=l,2,...;    T=1,...,m 
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and 


lim  Nj  =  2M  cos  9.  +  21^  cos  26 j ,    j=1,2, 
n->"» 


(6.48) 


Since  H^  =-  A  +  VN.V  (6. Ml),  one  can  write,  using  (6.6), 


lim  Hj  =  lim  (A  +  2B  cos  Q.    +  21^^  cos  29.)  ,   j  =  1,2,.. 


(6.49) 


This  means  that  lim^^„  H.,  j=l,2,...  have  the  same  configurations  as  A 
(6.3) 


v\rr\    \\.  -  \\ 


Tr\ 


n-»-< 


A-*' 


^<. 


^  ^Vi  ^ 


V\,j  Vn„^         \a,^ 


^i 

h.i 

^<i 

r 

V 

V>v> 

u., 

^» 

V 

v.> 

^o'i^r 

(6.50) 


m 


\  =  v^.--- 


vn 


*Note  limj^,„   ly^^    Nj    =  21 


m" 
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where 


h^.    =  p  +  2s  cos  e.  +  2  cos  26. 


(6.51) 


h^  .  =  q  +  2t  cos  9. 


LEMMA  6.4: 


lim   (Hj    -  UjAjUT)   =   Ojn    ,  j  =  l,2,...  (6.52) 


n-^oo 


Here 


(j)  (j) 

X^^   =   p  +   2s  cos   6.    +   2   cos  26::    +   2q  cos   6  +    4t   cos   6^    cos   6^        +   2r   cos  2 

j=1 ,2, . . .;    T=1 , . . .,m  (6.53) 


Finally,  we  obtain  the  following  result. 


THEOREM  6.1: 


lim  [a#  -  Z»A»Z,]  =  0(y'^)  (6.54) 

n-»-oo 


Here 


-5i|- 


Z^   =  X^d^xUj)  ,     Z.j  =  VU.jVUj  (6.55) 

m   in->-<» 


=  r  r  A.    r     i"^" 


Substituting  (6.4)  into  (6.53).  one  can  write 

(J) 
6. 


j^  =  16  [sin^  -2  +  Y^  sin2  _L_]  ,     t=1  ,2,  . .  .  ,m;  j=1.2,...   (6.56) 


6.  THEOREM  6.2;  Natural  frequencies  of  a  clamped  rectangular 
plate  have  the  following  explicit  expression 


"iT  =  •^  [ 5^^  i—W-'  J,T=1,2,...     (6.57) 

J^         1,2        n2     ''  nio 


Here  D  and  hIq  are  the  rigidity  of  a  thin  plate  and  mass  per  unit  area, 
respectively, 


e.  =^  e.  -  j  ,      e     =^e^-   ,    ,  (6.58) 


(t) 
9.  and  6^  are  zeros  of  A^   (9)   and  Aj^(9),   respectively.   They  are 

computed   for   large   m,n  >_  100.   Correction  terms  G.,e   satisfy 

inequalities 


0  <  £j,e^  <  .5  (6.59) 
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In  fact, 


I..   =  lim  _L_x|/2/- 


(j) 

r, .    ^         ■  2  ®j   . .    ^    .  2  ^T  1  /d~ 

=  [lim  sin  —  +  lim sin^  — _-J  y  — 

n^»  (Ax)2      2   ^^^  (^y)2       2  m^ 


[lim  ^ +  lim  —1 ]  /  — 

n>"  (n+1)^(Ax)^   m^"  (m+1)2(Ay)^    ""o 


TT^  J    +  1_]  /^  ,     j.T=1.2. 

&2        £2      m^ 


Note  that  eq.   (6.57)  generalizes  the  Navier  formula  (e^  =  e^  =  0)  for 
a  simply  supported  rectangular  plate  (see  Table). 
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rTTTTTTTTTTt 


£^-o 


6t,=  o 


O   ^£;   Z  ,15 


£^  =  o 
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TTTTTTTTTTTT 
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